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ABSTRACT 

Electrostatic Vlasov turbulence in a bounded spatial region is con- 
sidered. An iterative approximation method with a proof of convergence is 
constructed. The method is non-linear and applicable to strong turbulence. 

I. INTRODUCTION 

Consider the one dimensional Vlasov-Maxwell system of equations, 



— e»o 


for the electron distribution function, F(x,v,t), and the electric field, 
E(x,t), with a stationary and uniform ion background. Assune the existence 
of a solution of this system of equations for -1 i x i 1 (x is dimension- 
less, and measured in units of an arbitrarily length scale, L) , -® < v < •, 
and t > 0. Then a method for constructing approximations to this solution 
can be developed as follows: 
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The Fourier-Fourier transform method 1 can be used to transform 
equations 1-3 into an infinite system of first order hyperbolic partial 
differential equations. This system has been studied before. 1 One further 
aspect of this system will be considered here; it is that, on truncating 
the infinite system, the resulting finite system is of a standard form 
which has been used to produce constructive proofs of the existence of 
solutions to a wide class of such systems. 1 Here, the existence of a 
solution to the finite system will be assisted. The methods used for the 
proofs of existence will be applied, nevertheless, to produce approxi- 
mations to the exact solution of the finite system. The existence of the 
finite system solution and the issue of how well it approximates the 
solution of the infinite system will be addressed elsewhere. 

II. THE FOURIER-FOURIER TRANSFORM 
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then, f = F and e = E for -1 < x < 1, and on the boundaries 
f (±1, v ,t) = (1/2 [F( 1 ,v ,t ) + F(- 1 ,v ,t ) ] and e (±1 ,x ) = 

(1/2) [E(1,t ) + E(-1,t )]. Both f and e are periodic in x with period, 2. 
From equations 1 through 3, 




r 

« y 7 i- ‘ TTr 2 £„ \ 


H - -JO 


7. imTe* -- ^ - S J m (o ) Y) +■ |f M (, 1 - 


L_ ^ 

T > V 


r= o 


in which J_ = (-1) J , and 
m o 


where , 


XU, 7 ). ^ 

FU't.Vr- (], t‘" F U' ) v , ’ r 'l 


4 



Equations 7 and 8 are redundant when m i 0, and equation 7 yields no in- 
formation when m s 0. In the following equation 7 will be used to 
determine the e 0 ( T ) when m i 0, and equation 8 will be used to determine 

e 0 ( T >. 

It will be assumed in the following that f s 0 for m > M where M is 

n 

arbitrarily large, but finite. Since it is expected that those modes which 
have wavelengths comparable or shorter than the Debye length in the plasma 
will be strongly damped, 1 there is perhaps some a priori justification for 
the truncation. 

In Appendix A it is argued that the solution of equations 1-3 is deter- 
mined by a choice of F(x,v,0) on the initial plane (-1 ^ x 1, 

-• < v < ®, t = 0), of t ( v , t ) = F(-1,v,t ) - F(1,v,t ) in terms of F on 
the boundaries (x = ±.1, - • < v < ®, t > 0), and of c Q (0). That argunent 

does not depend on the assumption of a neutral plasma (space-averaged). In 
the following a neutral plasma will be assuned and then it will be shown 
that the choice of the equivalent quantities, f ffl (v,0), J q (v,t) and e Q (0) 
uniquely determines the solution of the truncated equations 6-8. 


III. THE NEUTRAL PLASMA 


The space- averaged electron density is f Q (0,T). The restriction to a 
neutral plasma is affected by setting f Q ( 0 , t ) s 1. This restriction is 
consistent with equation 6 only for a limited class of J q (v,t). From 
equation 6, 


HAM 

n 



where 6u(t) = u(-1,t ) - u(1,t ). In the following it will be assumed that 
J o (0,t) = 0, and that f Q (0,T,) = f Q (0,0) = 1. 
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Fyfe and Montgomery" have noted that e 0 (-r) cannot be chosen freely. 

They have produced an exact solution for e q (t) and u q (t) (the space average 
of u(x,t)) from their model of the one dimensional Vlasov-Maxwell plasma. 
Their results apply to the periodic plasma U(v,t) =0). A generalization 
to the non-periodic neutral plasma being considered here is possible. 

From equation 8 (or equation 3), 
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and from equation 6 (or equation 1), 
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where 6 P(t) = P (-1,t) - P (1,t) and 
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(Notice that since, for the neutral plasma, u(-1,t) = u(1,t), 6P ia 
actually just the difference in electron plasma temperature at the 
boundaries.) An exact solution of equations 9 and 10 is available; it is 
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The result of Fyfe and Montgomery is regained when 6P s 0. 

Notice that it is possible to obtain large e 0 (t) and u q (t) due to an 

approximately linear growth of the integral in equation 11 with increasing 

t if 6 P(t) contains harmonic oscillations with period 1 (the inverse plasma 

frequency). It is not possible for c q (t) to be constant in time unless 6P 

is also constant ,e Q (0) s 1/2 6P and u Q (0) s 0. Under these conditions 

u (t) 8 u rt (0) = 0. Since e (r) is a measure of the potential difference on 
0 0 o 

the boundaries, it should be noted that the preceding statements concerning 
c q (t) apply also to that potential difference. All of the above, and any 
other consequence of equation 11, apply exactly for the neutral plasma no 
matter what else is occur ing in the plasma. 


IV. BASIC. INXESRAl - EmJIQK 

Given the solution for e (t) and u (t), a major reduction in the com- 
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plexity of equations 6-8 can be made by introducing a new dependent 
variable through, y 
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then, 
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and the prime on the summation symbol indicates that the n = 0 term is 

omitted. Since e (t) can be considered a known function of time, a 
o 

solution of equation 12 for K m is equivalent to a solution of equations 6-8 
for f m> Notice that in the special case of a periodic plasma (J m = 0) 
equation 12 becomes independent of e Q (T). Thus, a single solution of 
equation 12, which will be shown to be determined solely by K ffl (v,0), is 
equivalent to the entire class of solutions for f which contains all 
possible choices of c q (0). 

Using the method of characteristics’ , equation 12 can be integrated 
once to obtain, 
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Equation 13 will play a central role In the following development of 


approximations to K . 

m 

V. THE APPROXIMATION METHOD 

The result which will be obtained in this section can be simply stated 
as follows: 

A sequence of functions, K (v,t;o), will be introduced with 

ID 
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Since F(x,v,0) will be assumed given, K m (v,0) can be considered a known 
function of v which is uniquely related to the initial F; a can be 
determined from a(v,t). The other members of the sequence are to be 
related to each other through, 
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It will be shown that K m ( v . t ; o) * K m (v,t); i.e., it will be shown that 
the sequence of approximations must converge to the exact solution. This 

convergence will not depend on the presence of ? small parameter for ex- 
pansion purposes and will apply for any finite value of M. Thus, this 
approximation method applies to strong turbulence with any finite number of 
wave modes, no matter how large. Notice that the character of the method 
is to place any member of the sequence in quadrature (through equation 16); 
it does not produce equations which must be solved. 

1 . Prelimina ries to Proof of Con vergence 

To facilitate the proof of convergence the K ffl will be asssumed vector 
components of a (2M + 1) - dimensional vector function, 

K(v,t) s (K_ m (v,t), K_ m (v,t) . . .K m (v, t) ) . Equation 16 can be considered 
an integral transformation which relates members of the sequence through 
K(o + 1) = T K(a ) . 

The proof of convergence will be given on the closed domain, 0, 
pictured in figure 1. An examination of equation 16 will show that 
knowledge of K(a) on 0 is necessary and sufficient to determine K(a ♦ 1 ) on 
0. The domain, 0, is centered on the line v s 0 since on that line all of 
the quantities of physical interest (the various moments of F as well as 
the Fourier components of the electric field) can be found. The time, T, 
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is ai.y finite time of interest; :he solution will be obtained for all 
0 < t i I, 

A further subdivision of 0 is necessary. Imagine the strips on the 
(v,t) - plane defined by r« < t (r ♦ 1) i where r s 0,1,2... and < > 0 is 
to be determined. Then, let 0^ be the intersection of the r' th strip with 
0. The transformation, T, will be shown to be a contraction* on each of 
the 0 p4 and then the results for the stationary element of T on each strip 
will be pieced together to yield equation 16. 

The following definition of a norm will be used. Let, 
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At each point, (v,t), the absolute values of all of the K 0 (v,t) are to be 
taken, and then the maximum of these chosen. Then, the supremum, on 0, of 
the resulting function is to be found. The vector, K, and the domain, 0, 
have been used here for illustrative purposes. Other vectors and domains 
will appear, but in each case, the symbol, || || Q , has the analogous 
meaning. 

In Appendix B it is shown that, in the limit M s •, ||K|| 0 = 1. In the 
following it will be assuned that ||K||q exists for finite M. This is not 
actually a new assumption; on 0, it is totally equivalent to the earlier 
assumption of the existence of a solution to the truncated system. It will 
also be assumed that the boundary conditions are chosen so that 
| |o| | Q exists. The number, N = 2(| |K| |q ♦ | | o I | Q T) , will be used. 
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2. Convergence, on a Narrow Strip 


Equation 13 can be used to show that 
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-ill now be given with K m (v,x; o) determined by equation 14. The symbol, 

T , will be used for the integral transformation in equation 17 on 0 . . 
re r 5 

The proof is in three parts . 
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3. Caigfirganfi&jai. large 


The proof of convergence on a narrow strip given in the preceding 
section will be used here to construct a proof of convergence at large (on 
the domain , 0) . It will be shown that , on °» iia K rn (v ' T; 0) = K m (V ’ T) 
where K m (v,x; o) is given by equation 14 and K o (v,x; a) is given by 
equation 16. An inductive argument will be given which assumes equation 18 
on 0 rfi as a starting point . 

Proof : 


Notice from equation 18 that, when r s 0, 
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for 0 < t i i, (It will be assumed that the values of v under consider- 
ation here are always on 0.) Assune, for some value of r, that, 
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b). Part Two 


Let V(v,t) and W(v,t) be any (2M + 1 )-dimensional vector fields on 0 rfi 

such that both 1 1 V|L i N and ||W|L IN. Define V' = T . V and 

u . u . ro 

r$ ro 

W' = T rfi W. Then, | |V» - W'| ! 0 1 1/2 | |V - W| l 0 for « small enough . 

r6 r6 

Proofs 
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Thu | | V * -W * | | Q i 1/2 ||V-W|| Q if 6 is small enough, 
rfi r« 

Plgcugsion 

This part of the proof shows that I . is a contraction. When T , is 

re ro 

applied to the difference of two normed vector fields, the resulting 

difference is reduced. In view of part a), of this proof, T r{ can be 

applied an arbitrary number of times to a pair of suitably chosen initial 

vector fields with the difference between the resulting vector fields 

reduced each time. In the following this basic property will be used to 

find the stationary element of T .. 

ro 

c). Part Three 

Let K m (v,T; o) be defined by equation 14. Define K m (v,T; a) s 

T„. K (v , t ; o-l) on 0„ .. Then, K (v,t) = lim K ( v,t; a) on 0„.. 
ro in ro id a+® in ro 

Proof 
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Therefore, ||K(0)|| o i (1/2) N <N. From part a), then, ||K(a)||_ N 

r6 r6 

for all a. Now, in the result of part b) . of this proof, let K(a) = V and 

K(B) = W. Then, 

lKW-KWl - 

' ' o rS o rj 

Suppose, a = B. Then, ||K(a) - K(B)||q = 0 for all a. Suppose, a > B 2 

r6 

1. Then, | |K(a) - K(B)| | Q l(1/2) 6 | |K(a -6) - K(0) | | Q i4N(1/2) 6 . 

rfi r6 

Thus, | |K(a) - K(B)| 0 as (a,B)+ ®. Similarly | |K(o) - K(B)||q-*" 0 as 

rfi r« 

(a, 8)- ® when B > a 2 1. Thus, K(a) is a Cauchy sequence 7 on 0 pfi . 

K(a) converges uniformly to K on 0 p5 where K* = K(a) 

• t « 

= lim T . K(a-1) = T .K . Since K = T . K it is a stationary element of 
a -►a* r ^ r o r 6 

T p ,. But, it is easy to see that the stationary element of T p g is unique. 

« i 

Suppose there are two stationary elements, K and L . Then, 
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Thus, | | K* - L*| | n = 0. 

r« 

Since K(v,t) is a stationary element of T p5 on 0 r{ , it is the station- 
ary element given by K(v,t) = lim K(v,t; a). Thus, equation 18 follows. 
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3. Convergence at large 


The proof of convergence on a narrow strip given in the preceding 
section will be used here to construct a proof of convergence at large (on 
the domain, 0). It will be shown that, on °» m K m (v ’ T * a) = K m (V ’ T) 
where K m (v,T; o) is given by equation 14 and K m (v,T; a) is given by 
equation 16, An inductive argument will be given which assumes equation 18 
on 0 pfi as a starting point. 

Proof; 


Notice from equation 18 that, when r : 0, 
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for 0 < t 1 {, (It will be assumed that the values of v under consider- 
ation here are always on 0.) Assune, for some value of r, that, 
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Notice from equation 19 that equation 20 is true when r = 1. By the 
induction hypothesis (equation 20), 

I* 

0 



Equation 21 can be substituted into equation 18 to obtain, 
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on 0 r5 ; i.e., for rfi < t < (r + 1) 6. In particular, equation 22 is true 
for t = (r + 1) 6. Thus, if equation 20 is true for any value of r, it is 
true for all larger values of r. Since equation 20 is true for r = 1, it 
is true for all values of r and equation 22 is true everywhere on 0. 
Discussion : 

Equation 22 is the primary result of this paper. From equation 22 a 
sequence of functions can be computed with the understanding that the 
sequence will converge to the truncated Fourier series expansion of the 
Vlasov plasma distribution function. 
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VI. CONVERGING SEQUENCES FOR. THE. KINETIC AND FIELD ENERGIES 


In Appendix B it is shown that, in the limit M = », K' and K" (K* = 

m mm 

3K ffl /3v, etc.) can be uniformly bounded on 0. In the following it will be 
assumed tha* | I K * | | Q amd ||K"||q exist for finite M. In this case essen- 
tially the same procedure as given above can be carried out to prove, 
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and, 
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as long as the boundary conditions are chosen so that | | a * | | Q and 
I I a " | | 0 exist. 
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Equations 23 and 24 show that the sequences of functions which are 

obtained by differentiating the K m (v,T; «) converge to the respective 

derivatives of K (v,t). Thus, K' (a) and K" (a) can be expected to approx- 
m m m 

imate the exact derivatives. The derivatives, with respect to v, play an 
important role in the application of this theory. At v = 0, K ffl and its 
derivatives are related to the coefficients in Fourier series expansions of 
the moments of the distribution function. For example, 


tx> 
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do 
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in which f and f’ can be computed from K m and K* (see Appendix B) . 
From equations 23 and 24 it can be shown that, 
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This equation bears on energy conservation in the plasma. Using 
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and P Q ( t ) for the space average of P(x,t), equation 25 can be rewritten in 
terms of more familiar notation as, 

M / 

P 0 U \*\ z-t W + HA* )_ fef U n t»,e)| 
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where 6Q(x) = Q( — 1 , X ) - Q(1,X). The a'th iterate to the space averaged 

kinetic energy density is given by P o (t; a). The contributions to 

P q (t;o) on the right side of equation 26 are as follows: The first line in 

equation 26 gives the total energy in the plasma at t = 0. P Q (0) is the 

2 

initial kinetic energy, c Q (0) is the initial electric field energy in the 
space averaged part of the field, and the sun on this line represents the 
field energy in E'(x,0) averaged over space. The second line of equation 
26 gives the negative of the total field energy at any time. The sun in 
the second line is the space average of the field energy in the turbulent 
part (E'(x,t)) of the electric field, as given by the (o-l)’th iterate. 

The last line of equation 26 gives the accunulated net transfer of energy 
into -1 < x < 1. Thus, energy is conserved at each iteration. Further- 
more, sequences of approximations for the kinetic and electric field 
energies in the plasma can be computed from equation 26. From equations 23 
and 24 it can be seen that these sequences must converge to, respectively, 
the kinetic energy in the truncated distribution function and the 
corresponding turbulent electric field energy. 
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VII. CONCLUSION 


A one dimensional electrostatic Vlasov-Maxwell plasma model has been 
considered in a bounded spatial region. Consideration has been limited to 
a plasma with uniform and stationary ion background and with zero 
space-averaged charge. An iterative method has beon constructed for 
computing a sequence of approximations to the probability distribution 
function for the initial-boundary value problem. 

The probability distribution function has been Fourier transformed in 
its velocity variable, and Fourier series expanded in its spatial variable. 
The Fourier series expansion has been truncated at an arbitrarily large but 
finite value. It has been assisted that a solution exists to the finite 
system of partial differential equations which govern the truncated ex- 
pansion of the distribution function. Under this assumption it has been 
shown that the sequence of approximations mentioned in the preceding 
paragraph must converge to the exact solution of the truncated system. 
Convergence does not depend on the presence of a small expansion parameter 
for expansion purposes. The degree to which the solution of the truncated 
system approximates the solution of the infinite system has not been 
considered, but, in view of the arbitrarily large number of Fourier modes 
that is allowed, a good approximation is anticipated in many applications. 

The issue of the rate at which convergence occurs is under investi- 
gation at present. In those situations where convergence is rapid enough 
to make this iterative method useful, it can be viewed as an approximation 
technique for Vlasov turbulence which is non-linear and applicable to 
strong turbulence in a bounded region of space. 
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AEJEEMfiULA 


Consider the one-dimensional Vlasov-Msxwell system of equations 1-3. 
Assune the existence of a solution to this system of equations on the 
domain, D, defined by -1 < x < 1, - • < v < •, and t > 0. Then, what 
combination of initial and boundary conditions on F and E uniquely, 
consistently, and conveniently determines that solution? 

Given the existence of E(x,t), the method of characteristics can be 
used to solve equation 1.* The solution is F(x(s), v(s), t(s)) * 

F(x(0), v(0), t ( 0) ) where x(s), v(s), and t(s) are the solutions of the 
system, 


cUl*) 

-r - v 

cU 

iUil' - - 

(A ^ 

Jl’iw , 

TT = 1 

subject to x(0), v(0), and t( 0) for initial conditions. In typical 
applications of this type of solution, the point, (x(0) , v(0) , t( 0)), is 
assumed on one of the boundaries where F(x(0), v(0), t(0)) is known,. 
Then, F is known everywhere along the characteristic line given by x(s), 
v(s), and t(s) for s 2. 0. The complete solution depends on filling all 
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(x,v,x) of interest with characteristic lines which are connected to 
boundaries where F(x(0), v(0),x(0)) is known. But, this typical approach 
is self-contradictory for the initial-boundary value problem being 
considered here. 

Figure 2 contains schematic representations of projections onto the 
plane, v s 0, of various types of possible characteristic lines. 

Progression along the characteristic lines, with s increasing, is indicated 
by arrows. The dashed lines indicate characteristics which enter D on the 
initial plane (-1 S. x < 1, -• < v < ®, t = 0) or on either of the 
boundaries (x = ±1 , x >0) and then remain trapped in D. 

These are the only characteristics that can be treated as outlined in the 
preceding paragraph. All other characteristic line*- (solid lines) enter D 
and then exit , with s increasing. The solution of equation 1 gives F = 
constant along each of these characteristic lines. Thus, F must have a 
single value at every pair of entry and exit points. The initial and 
boundary values for F are not independent and cannot be chosen arbitrarily. 

There is no unique method for choosing the initial -bound ary values for 
F such that the possible contradictions discussed above are avoided. The 
method used here has been chosen for both mathematical and observational 
convenience. In figure 3 it is demonstrated that the solution at an 
arbitrary point, (x,v,x) , is determined by specifying F(x,v,0) on the 
initial plane, and a(v,t) s F(-1,v,t) - F(1,v,x) in terms of F on the 
boundaries. The solid lines indicate possible characteristic lines along 
which F = constant . The dashed lines connect exit and entry points at 
which the value of F must be related through the use of A(v,x). In the 
example presented in figure 3, F(x,v,x) = F(x Q . v Q , 0) 

+ a(v.j, t^) - A(Vg> tg). Given the existence of E(x,x), this choice of 
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F(x,v,0) and a(v,t) uniquely determines F(x,v,t) through equation 1. The 
type of contradictions discussed above are avoided because A(v,t) always 
relates F on an exiting characteristic line to F on an entering 
characteristic line. 

The A-method for specifying the boundary data allows for the expression 
of the data in terms of the measurable moments of the distribution function 
on the boundaries (e.g., the difference in density, current, temperature, 
etc.,) where others do not. (Mathematically , it is acceptable to specify F 
on the boundaries for incoming velocities only. It is also possible to 
express the boundary conditions in terms of the ratio, F(x * 1)/F(x * -1). 
Neither of these boundary conditions can be expressed in terms of the 
moments on tne distribution function on the boundaries.) Notice that the 
well studied "periodic plasma" in which A(v,t) : 0 is a special case of the 
A -method. More generally, the A -method plays a natural role in the Fourier 
series expansion analysis of the non-periodic plasma. A choice of F(x,v, 

0) and A(v,t) leads to a unique set of F m (v,0) ard J m (v ,t ) . Thus, the 
A -method has been chosen as the basis for this study of the 
initial-boundary value problem for the Vlasov-Maxwell plasma. 

Now, assume F(x,v,t) is known. Then, what additional initial-boundary 
data must be specified to determine E(x,t) from equations 2 and 3? 

If E(x,t) is separated into its space average, c Q (t), plus an 

I 

x-dependent part through E(x,t) * e o (t) * E (x,t), then equation 2 

t 

determines E (x,t) only. From equation 2, 
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- 1 


where n Q (T) is the space average of n(x,T), and 


A. 2 


^ \\) - \ X 


From equation 3. 


A. 3 



U 0 rO 


where u q (t) is the space average of u(x,t). Notice, if F is known, then 
♦ 

E (x,t) is given by equation A .1 with no freedom for choosing boundary or 
initial conditions. Further, given F, u q (t) can be calculated, and then, 
c q (t) can be calculated from equation A. 3 if e Q (0) is specified. Thus, 
only c q (0) need be specified to determine E(x,t) given F. 

By combining the argument given above for determining F given E and for 
determining E given F, it seems plausible that specification of F(x,v,0), 
a(v,t) and e o (0) uniquely determines a solution of equations 1-3. This is 
not a proof of the existence, nor the uniqueness of that solution; it is, 
at best, a plausibility argunent . 
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APPENDIX B 


The Fourier series expansion of F has been truncated in section II, and 
the resulting truncated version of equations 6-8 has been solved in 
principle in section VI, not for F, but for 


B.l 


'T~ c *1 “ * i ( , - c >1 v' V r , 

-- /_ e 

m - - il 


— CKJ 


The f m (v,-r) sre the solution of the truncated equations which are still 
related to f(x,v,x) in the usual manner: 




B.2 


1 






CK> 


The degree to which f approximates F has not been investigated in this 
paper. However, since the theory dev*. 1 ped in this paper applies to a 
truncation for arbitrarily large M it has been assumed that this 
approximation can be made as good as necessary. In particular, it has been 
assumed that for some finite M, ||K|( 0 , ||K|| 0 and | | K" | | Q exist . The 
conditions under which this assumption is valid are under investigation. 

The goal of this appendix is to prove that ||K|| 0 = 1 and both 

t i 

|| K | | q and | |K | |q exist when M = ®. For this purpose it is convenient 
to introduce 
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B.3 
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In the following it will be assumed that the boundary conditions are chosen 
such that, when the solution of equation 11 is substituted into equation 
B.3, the resulting 4>(x) can be uniformly bounded on the interval, 0 < t 
T. From the definition of the K m , 


B -' 1 I'M" l VJ 

\K\ *- T 141 1U + I V m \ 


and 

u:|i ici 

» ft 

If f m’ f m and f m can be uniforml y bounded on 0, then the goal of this 
appendix will have been achieved. 

From equation B.2, 
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But, in section III it was shown that f Q (0,T) = 1, Thus |f (v,t)| <. 1 and 

I 1 K| | Q = 1. 

By differentiating equation B.2 it can be shown that, 



» 2 n 

where f q (0,t) < 0. But, - (1/ir) f q (0,t) is actually the kinetic energy 
in the plasma. The equation which governs conservation of energy can be 
obtained from equations 6-8 and written, 



n-- 


y 

where 6Q(X) is the net rate at which energy is entering the region in x of 


interest (see section VII). 



Thus , 

• (ft* c « to ' 1 + 



V|S -OO 
Ks 

+ i \llx l Q\\) 

The first line of the right side of this equation represents the total 
initial energy in the plasma. If this initial energy is chosen bounded, 
and if the rate at which energy is allowed to enter, 6 Q(t), is assumed 
uniformly bounded on 0 i t i T, then |f" n (v,T)| is uniformly bounded on 0. 


Given the uniform bounds of f o (0,T) and f" o (0,T) which have been 

obtained above, it is, possible to prove a uniform bound on f* . From 

m 

equation B.2, ^ 

I CM - * 

- » Jaa 

* i tfo 

- "a \^\ \‘^ v * 

— \ - co -i » 



1 -If 


Thus, f is uniformly bounded on 0, as well as f and f , and from 

ID IQ ID ' 

f ft 

equations B.5 and B.6 | |K | | Q and | |K | | Q must exist. 
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FIGURE 1 


The domain 0 on the (v,t) -plane. 


FIGURE 2 A schematic representation of all possible characteristic 

curves on the domain D. 

FIGURE 3 A possible construction of F(x,v,t) using F(x,v,0) and 

a(v,t) . 
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